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This paper, which is a continuation of two earlier papers of the same title, contains 
numerical results for the field anomaly near a coastline when the surface impedance changes 
in a linear manner between land and sea. The earlier results for an abrupt boundary are 
recovered as the width of the transition region is reduced to zero. In general, it is found 
that the characteristics of the transition region will not produce significant modifications of 
the transmitted field. However, the magnitude of the reflected field is greatly reduced as 
the width of the transition zone is increased beyond about one-quarter wavelength. 



1. Introduction 

In part I [Wait, 1963] of a series of papers of this 
title, the propagation of radio waves across a flat- 
lying coastline was investigated theoretically. In 
part II [Wait and Jackson, 1963], the influence of an 
elevation change between land and sea was con- 
sidered. It is the purpose of the present paper to 
consider again the flat-lying coastline but with special 
attention given to a nonabrupt variation of the 
conductivity at the coastline. 

As indicated in part I, the assumption of a sudden 
change of electrical properties gives rise to an appar- 
ent singularity of the field right at the coastline. It 
is shown in what follows that the " singularity' ' is not 
present when the surface impedance changes grad- 
ually at the junction or coastline. 

2. Formulation 

The model is the same as that used in part I. The 
essential features are described briefly here. The 
transmitting antenna, A, located on the land, is at a 
great distance from the coastline. For present pur- 
poses, it is only necessary to assume that A is a source 
of vertically polarized groundwaves. The receiving 
antenna, B, is relatively near the coastline but it may 
be on the sea. It is further assumed that B is equiva- 
lent to an infinitesimal vertical dipole and thus it 
responds only to the vertical electric field. 

From consideration of reciprocity, it is clear that 
the role of transmitter and receiver may be inter- 
changed. Thus, in general, it is most meaningful to 
express the results in terms of mutual impedance 
z m -\-Az m between the respective terminal pairs of 
antennas A and B. Here z m is the mutual impedance 
if the surface impedance were a constant Z jov all 
points of the earth's surface. Thus, Az m is the 
modification of the mutual impedance which results 
from the presence of the inhomogeneity of surface 
impedance Z' '. In the present problem, the surface 



impedance of the land is Z while that of the sea is 
Z\ The situation is illustrated in figure 1 where the 
earth's surface is the (xy) plane of a Cartesian coor- 
dinate system and the coastline is at x=0. 

3. Statement of Formulas 

Quoting from part I, the working formula for the 
mutual impedance increment Az m is given by 



Az m % ,„ , 
— ^— ~ e lkC i d i 

Zm, £ 



where 



£[c/-f)-^ ( x)],--, 



Xm 2) [kC\\x-<h\}dx, (1) 



f(x) = (Z f -Z)/ Vo , ?7o=12(K 

H ^ is the Hankel function of order zero of the 

second kind, 
Ci = cos d , &=27r/wavelength, 
di = shortest distance from coastline to B (positive if 

B is over sea and negative if B is overland), 

and, finally, e is chosen to be sufficiently large that 
f(x)^0 for arC— e. 

For purposes of the present paper, some further 
simplifications may be made. Specifically, it is 
noted that 

f(x)=-Ao(x)e i "\ (2) 

where A (x) is approximately a real positive quantity 
whose magnitude is small compared with unity (e.g., 
A <^0.1). For example, at sufficiently large positive 
values of x, 

A (x)^A ^( — j-. ) ( 1 XT — ) 

( €o0) V7 iepa \ 1/2 ( 
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Figure la. Plan view of the earth's surface] illustrating) a 
coastline with a finite transition zone between land and sea. 
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Figure lb, c. The linear surface impedance variation. 

where <r g and e g are the conductivity and permittivity 
of the land and a and e are the conductivity and 
permittivity of the sea. In most cases of practical 
interest, co is sufficiently low and a' is sufficiently 
large that 

(4) 



Ao^(W*«) 1/2 , 



which is real. 



4. Linear Variation of Surface Impedance 

In part I, it was assumed that A (x) was a step 
function at x=0 corresponding to an abrupt bound- 
ary. Here it is desirable to allow A (x) to have a 



finite transition region of width d where A (x) varies 
in a linear manner. Thus 

Ao(aO=Oforx<—do/2, 

+^]A for-r/ /2<x<r/ /2 ) 
=A for x>d /2. (5) 



-p 



With this model, the surface impedance is imagined 
to change in a linear manner over the distance d as 
indicated in figure lb. 

To simplify the computational problem, certain 
dimensionless quantities are introduced as follows. 

a=W 1 x, a 1 =kC 1 d u D=kC 1 d , 

and D =D/C 1 =kd . Then, it easily follows that (1) 
may be written in the form 



Az m A e* 3W4 



£ 



H^[\a- ai \]e-^ 



(6) 



where, in the case of a linear variation, 
G=0iora<-D/2, 

_ a +(D/2) 



D 
-I fora>D/2. 



for -£>/2<a<.D/2, 



(7) 



Therefore, 



^=i for -D/2<a<D/2, 
da U 



= for a<-D/2 and a>#/2, 



(8) 



as indicated in figure lc where G is plotted versus a. 
The mutual impedance formula now becomes 



%m L ( lJ 



where 



Qi- 



9 z37r/4 



2D 



J -D, 



HP\\a-a, 



(9) 



(10) 



and 



,ix/4 r rD/2 a A-T)\2 



■ai\\e' 



ia do 



J D/2 



(2)1 



a—oLi\\e 



e~ ia da V- 



(11) 
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To effecl these integrations, il is desirable to make 
use of the following two identities: 



d_ 

il.r 



[xe± tz [m 2) (x)^iH™(x)]]=e* 



'WW, (12) 



and 



^Tj^[xH?>(x) + (lTix) 



f\x)]\=xe± iz m n {x), 
(13) 



where upper (or lower) signs are to be considered 
together. In verifying these identities, it is neces- 
sary to employ two well-known relations in Bessel 
function theory [MacLachlan, 1934]. These are 



d_ 

dx 



m*>(x)=-m*>(x), 



and 



if H {2) ( r) 

-f m 2) (x)=m 2) (x)-thA£i 

(IX x 



(14) 



(15) 



It is now apparent that Azjz m may be expressed 
entirely in terms of Hankel functions of order zero 
and one, with various real arguments. Using such 
expressions, curves of the real and imaginary parts 
of Q have been prepared. Some of these are shown 
in figures 2a to Sb where, in each case, the abscissa 
is the parameter a x or kCidi. 

For purposes of discussion, it is convenient to de- 
scribe Q as the field anomaly resulting from the 
inhomogeneity. In view of the relation 



^=A [Re Q+i Im Q], 



(16) 



it is evident, for A real, that Re Q is a measure of the 
amplitude change of the field where Im Q is a meas- 
ure of the phase. While the results are still valid if 
A is complex, the above simple interpretation is not 
applicable. Also, it should be remembered that under 
all conditions, |As OT /2 TO |«l or |A (?|<<1, if the pres- 
ent results are to be given any confidence. 

5. Discussion of Numerical Results 

The curves in figures 2a and 2b are applicable to 
normal incidence. The sinusoidal-like ripples for 
negative values of a x may be regarded as an interfer- 
ence pattern resulting from the combination of the 
incident wave and the reflected wave. It is apparent 
that, as the electrical width D of the transition zone 
is increased, the magnitude of the reflected wave is 
generally decreased. On the other hand, the char- 
acteristics of the transmitted wave, for large positive 
values of <x u are not appreciably modified. However, 
it may be observed that in the proximity of the coast- 
line, the nature of the field is profoundly influenced by 
the width of the transition region. In general, the 
rapid variations of the field are smoothed out when 
the transition distance D is increased. As may be 




Figure 2. The real and the imaginary parts of Q as a function 
of <x\ or KCi (\\ where dj is the distance from the coastline. 

seen in figure 1, the field is beginning to have a singu- 
lar behavior for the smallest value of D Q (i.e., 1) which 
is very similar to the corresponding curves given in 
part I for the abrupt boundary (i.e., D =0). 

The curves in figures 3a and 3b for 6 =20° are very 
similar to those in figures 2a and 2b for d =0°. Evi- 
dently, a slight obliquity in the traverse across the 
coastline does not modify noticeably the field be- 
havior. However, for more oblique conditions, the 
curves in figures 4a, 4b, 5a, and 5b for O =4O° and 50° 
indicate that significant modifications take place for 
angles in this range. Principally, it is noted that 
when #o is in the vicinity of 45° the reflected wave is 
greatly reduced in magnitude. This effect was also 
present in the case of an abrupt boundary. 

The curves in figures 5a to Sb illustrate the be- 
havior of the field anomaly for rather highly oblique 
angles. Again, as in the case of the abrupt boundary, 
the reflected wave is quite strong. Furthermore, I he 
magnitude of the reflected wave is not appreciably 
reduced when D is increased from 1 to 5. However, 
in nearly all angles O the behavior of the transmitted 
field is not appreciably modified by changes in the 
width of the transition zone. 
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Figure 3. The real and the imaginary parts of Q as a function 
of a\ or KCi dj where dj is the distance from the coastline. 



The relative insensitivity of the Q versus a x curves, 
at highly oblique angles, to variations in D is reminis- 
cent of the reflection of waves from horizontally 
stratified media [Wait, 1962]. In the latter case it is 
known that the influence of diffusiveness or non- 
sharpness of the boundary is minimized at highly 
oblique incidence. 

6, Conclusions 

The results given in this paper would seem to shed 
considerable light on the nature of the vertical elec- 
tric field near the boundary of separation. The 
singularly or infinite behavior of the field obtained 
in earlier studies, is not present when the surface 
impedance changes gradually between the land and 
sea portions. However, as the transition zone is 
diminished the field becomes rapidly varying in the 
vicinity of the coastline and becomes very similar to 
the predicted behavior for an abrupt boundary. A 
conclusion similar to this was arrived at by Godzinski 
[1962] whose analysis communicated to me in outline 
was restricted to normal incidence (i.e., O =O°). 




Figure 4. The real and the imaginary parts of Q as a function 
of a\ or KCq di where dj is the distance from the coastline 




Figure 5. The real and the imaginary parts of Q as a function 
of a\ or KCi di where d\ is the distance from the coastline. 
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Figure 6. The real and the imaginary parts of Q as a function 
of a\ or KC^ di where d\ is the distance from the coastline. 
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i The careful drafting of the figures is due to Jerry Hodges. 
2vz/c should be replaced byc/tfxineqs (29) and (31), v y /c should be replaced by 
c/Vy in eqs (30) and (32), and v y jv x should be replaced by v x /v y in eq (33). 
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Figure 7. The real and the imaginary parts of Q as a function 
of <x\ or KCj di where d\ is the distance from the coastline. 
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Figure 8. The real and the imaginary parts of Q as a function 
of a\ or KCi di where dj is the distance from the coastline. 
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